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1 Introduction 



The phase diagram of QCD in the density (pg) - temperature (T) plane has been ex- 
plored by many authors; QGP in high-T region or color superconductivity in high-p^ 
region is a typical phase in that plane. Here we are interested in low temperature and 
moderate density region relevant to compact stars, where magnetic order is expected. 

Origin of the magnetic field in compact stars is one of the long-standing problems 
since the first discovery of pulsars in early seventies. Recent discovery of magne- 
tars with huge magnetic field of O(10 14_15 G)has revived this problem. Since nuclear 
matter is developed inside compact stars, we are tempted to consider spontaneous 
spin polarization of nucleons as a microscopic origin of the magnetic field. Realistic 
calculations have been performd for polarized nuclear matter, but they lead us to the 
negative results pQ. In ref.[2] possibility of spin polarization of quark matter has been 
suggested by a simple consideration in analogy with itinerant electrons, where the 
Fock exchange interaction is responsible to ferromagnetism [3]. A weakly first-order 
phase transition has been demonstrated around nuclear density. Using this result 
we can roughly estimate the magnitude of the magnetic field at the surface to be 
O(10 15 ~ 17 G), which may explain the magnetic field of magnetars. The coexistence of 
ferromagnetism with color superconductivity has been also discussed in ref. [3]. 

Here we apply the Landau Fermi-liquid theory (FLT) to elucidate the critical 
behavior of the magnetic phase transition at finite density and temperature [5]. We 
evaluate the magnetic susceptibility of quark matter. The divergence and sign change 
of the magnetic susceptibility is a signal of the magnetic instability to the ferromag- 
netic phase, since its inverse measures the curvature of the free energy at the origin 
with respect to the magnetization. Thus quarks near the Fermi surface are respon- 
sible to the magnetic transition and the spin dependent quark-quark interaction and 
the density of states near the Fermi surface are the key ingredients within FLT. 
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Theoretically we find a non-Fermi-liquid behavior of the magnetic susceptibility.lt 
is well known that there appears a non-Fermi-liquid behavior in the expression of the 
specific heat in QCD as well as QED, which is caused by the transverse gauge field 
because it is not statically screened [6]. 

2 Relativistic Fermi liquid theory 

Within the Landau Fermi-liquid theory (FLT) we assume a one-to-one correspondence 
between the states of the free Fermi gas and those of the interacting system [5] . Quarks 
are treated as quasi-particles carrying the same quantum numbers of the free quarks, 
and the quasi-particle distribution function is simply given by the Fermi-Dirac one, 

n(k,C) = [l+exp(/3(e k , c - / i))]- 1 (1) 

with the quasi-particle energy specified by the momentum k and a spin quantum 
number ( = ±1. 

2.1 Screening effect on the quasi-particle interaction 

In the following we consider the color-symmetric interaction among quasi-particles 
that can be written as the sum of two parts, the spin independent (/k,q) an d dependent 
(/k, q ) terms; 

/kC,qC = /k,q + CC7k,q- ( 2 ) 

Since quark matter is color singlet as a whole, the Fock exchange interaction gives 
a leading contribution. We, hereafter, consider the one-gluon-exchange interaction 
(OGE). For a pair with color index (a,b), the Fock exchange interaction gives a 
factor (X a )ab(^a)ba = 1/2 ~ l/(2iV c )£ a fc, which is always positive for any pair. Hence 
the situation is very similar to electron gas in QED. Since we are interested in the 
electromagnetic properties of quark matter, only the color symmetric interaction is 
relevant, which is written as 

1 ^ m m 

c a,b ^ 1 

with the invariant matrix element, 

M kf , qC , = -9 2 ^tr (XJ2XJ2) M^(k, C; q, C)D^(k - q), (4) 

where M^(k, C; q, C) = tr [fp{k, ChV(<7, CO]- 

Since the OGE interaction is a long-range force and we consider the small energy- 
momentum transfer between quasi-particles, we must treat the gluon propagator by 



2 



taking into account HDL resummation. Thus we take into account the screening 
effect, 

D^k -Q)= KMp) + P^A(P) - ^ (5) 

with p — k — q, where D t ^)(p) = (p 2 — IL^)) -1 , and the last term represents the gauge 
dependence with a parameter £. Pfffi is the projection operator onto the transverse 
(longitudinal) mode, 

P% = (l-^o)(l-^o)(-^-^) 

PU = -^ + Vj f-?l- (6) 
The self-energies for the transverse and longitudinal gluons are given as 

n,,„ P ) - + 

n«(Po,P) = -i E (7) 



f=u,d,s 



PI 



in the limit po/|p| — *■ 0, with ufj = kpj/Epj and the Debye mass for each flavor, 
m|)j = g 2 fXfkF,f/2Tc 2 . Thus the longitudinal gluons are statically screened to have 
the Debye mass, while the transverse gluons are dynamically screened by the Landau 
damping, in the limit Po/|p| ~~ 0. Accordingly, the screening effect for the transverse 
gluons is ineffective at T = 0, where soft gluons (po/|p| ~~ 0) contribute. At finite 
temperature, gluons with p ~ 0(T) can contribute due to the diffuseness of the 
Fermi surface and the transverse gluons are effectively screened. 



2.2 Magnetic susceptibility 

We consider the linear response of the normal (unpolarized) quark matter by applying 
a small magnetic field B. Using the Gordon identity, the coupling term with the 
uniform magnetic field (A = B x r/2) can be written as 

J d A xC- mt = e q J d A xq~f ■ Aq 

= fi q Jd 4 xq[L + i:}-Bq, (8) 

with /j, q being the Dirac magneton. We discard the contribution of the orbital angular 
momentum (L) by assuming the uniform distribution of quarks. Thus the magneti- 
zation (M) / for each flavor can be written as 

(M) f = V- l (Jd 3 xq f i: z q f ), (9) 
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where we take B//z. Accordingly the magnetic susceptibility is defined as 



v _ V v' - V d{M) f 
XM - 2^ Xm~ 1^ a R 

j=u,d,s j=u,a,s 



(10) 

5=0 



Hereafter, we shall concentrate on one flavor and omit the flavor indices because the 
magnetic susceptibility is given by the sum of the contribution from each flavors. The 
magnetic susceptibility is proportional to the number difference between different spin 
states (C = ±1); it is explicitly caused by the applying magnetic field, and implicitly 
caused through the spin-dependent interaction, 

6n kf=+ i - 5n k c=-i = {-QD^qB + 5e kC=+1 - <5e kC= _i] (11) 

<7e k 

with the gyromagnetic ratio go ~ 2, where 

k( = ^c E / ^^/kCiqC'^qC- ( 12 ) 



5e 



(2 

Magnetic susceptibility is then written in terms of the quasi-particle interaction 



fgn^V N(T) 
Xm - [—) l + N{T)Ja (13) 

where gp> is an angle average of go, and f a is the Landau-Migdal parameter averaged 
over the Fermi surface [HJ [12] . 



3 Magnetic properties at T = 

N(T) is the effective density of states at the Fermi surface, and is simply written as 



N-\0) 



7T 



N c k 2 F 



(14) 



in the limit of zero temperature . Eq. (fT4|) defines the Fermi velocity, which is given 
by using the Lorentz transformation [5], 



v F 



<9n k 



kp 

|k|=fc F ^ 



N h 2 
9_2 '1' 



(15) 



where ff is a spin-averaged Landau-Migdal parameter. 

Finally the magnetic susceptibility at zero temperature can be written in terms 
of the Landau-Migdal parameters, 



Xm — XPauii 



1 + 



N c k F fx 



7T- 



\fi + r 



(16) 
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where XPauii is the usual one for the Pauli paramagnetism, XPauii = 9D^ q N c kF^/ ^ ■ 
The quasiparticle interaction on the Fermi surface can be written as 



qC'l|k| = |q|=fcj: 



-a 



m 

9 M 



-M m D L (k - q) + M u D T (k - q) 



(17) 



with the effective coupling strength, C g = N £ N 2 g 2 - 

We can see that the both Landau parameters /{, f a include the infrared singular- 
ities due to the absence of the static screening for the transverse gluons; Dx(k — q) ~ 
— l/(k — q) 2 = —1/2^(1 — cos^q) in this case, so that the logarithmic diver- 
gences appear in the Landau parameters through the integral over the relative angle, 
J^ kq l/(l-cos^ kq ). 

Finally magnetic susceptibility is given as a sum of the contributions of the bare 
interaction and the static screening effect. We can see that the logarithmic divergences 
exactly cancel each other to give a finite result for susceptibility [HI E]. 



[Xm/ XPauli) 



-i 



CgN c fl 



12n 2 Elk 



F^F 



m(2E I 



m) 



4E F m 



2m 2 )«;m — 

K 



with k = m 2 D /2k 2 F . Obviously this expression is reduced to the simple OGE case 
without screening in the limit k — > 0; one can see that the interaction among massless 
quarks gives a null contribution for the magnetic transition. The effect of the static 
screening for the longitudinal gluons gives the contribution of g A ln(l/g 2 ). In the 
nonrelativistic limit, it recovers the corresponding term in the RPA calculation of 
electron gas [31 El EJ - 
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Figure 1: Magnetic susceptibility at T = 0. The solid curve shows the result using 
simple OGE, while the dashed and dash-dotted ones show the screening effects with 
Nf = l(only s quarks) and with Nf = 3(u, d, and s quarks) respectively. 
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In Fig. 1, we plot the magnetic susceptibility at T = [8], [9j. We take the QCD 
coupling constant as a s = g 2 j^n = 2.2 and the strange quark mass m s = 300MeV in- 
ferred from the MIT bag model. We consider here the MIT bag model as an effective 
model succeeded in reproducing the low-lying hadron spectra. The coupling constant 
looks rather large, but this value is required for the color magnetic interaction to ex- 
plain the mass splitting of hadrons with different spins; e.g. for nucleon and A isobar. 
We think this feature is relevant in our study, because the coupling constant is closely 
related to the strength of the spin-spin interaction between quarks in this model. 
Moreover, the quark density in the MIT bag model is moderate, 0.25fm -3 , which is 
the similar one we are interested in. Note that the perturbation method should be 
still meaningful even for this rather large coupling, since the renormalization-group 
analysis has shown that the relevant expansion parameter is not the gauge coupling 
constant g 2 but the product of g 2 with the Fermi velocity vp, which always goes to 
zero as one approaches to the Fermi surface [7]. 

One can see that the magnetic susceptibility for the simple OGE without screening 
diverges around kp = 1.3 fm . This is consistent with the previous result for the 
energy calculation. [HIS]- One may expect that the screening effect weakens the Fock 
exchange interaction so that the critical density get lower once we take into account 
the screening effect. However, this is not necessarily the case in QCD. The screening 
effect behaves in different ways depending on the number of flavors. Compare the 
results for the Nf = 3 with the one for Nf = 1. In the case of Nf = 1, k < 2 the 
screening effect works against the magnetic phase transition as in QED. However, for 
Nf = 3, k > 2 so that the critical density is increased. Consequently the screening 
effect does not necessarily work against the magnetic instability, which is a different 
aspect from electron gas [9]. 



4 Finite temperature effects and Non- Fermi- liquid 
behavior 

At finite temperature, the magnetic susceptibility is given by 



where /" and /" denote the longitudinal and transverse parts of f a respectively 
[H1H2]- First, we evaluate the effective density of states on the Fermi surface defined 




(19) 



by 




(20) 
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with eo = e|k|=o- The quasi-particle energy u should be given as a solution of the 
equation, 

w = E ki u) + ReS+(u;, k(u)), (21) 

where we discard the imaginary part within the quasi-particle approximation. 

The one-loop self-energy is almost independent of the momentum, and can be 
written as HU1 



ReE+(u;, k) ~ ReE+(/x, k F ) - C {f U / (to - y) In , A , + A rcg (w - fj) (22) 

12tv z \uj — fj,\ 

around u ~ /x with C/ = (JV^ — l)/(2iV c ) and = k F /E kp . A is a cut-off factor 
and should be an order of the Debye mass, A ~ 0(m£>). Note that the anomalous 
term in Eq. ( 1221) appears from the dynamic screening of the transverse gluons, and 
the contribution by the longitudinal gluons is summarized in the regular function 
A reg (cj — /j) of 0(g 2 ). Within the approximation given by Eqs. (1221) and (??), the 
self-energy is independent of spatial momentum k and thus we omit the argument k 
hereafter. The renormalization factor z + (k) is then given by the equation, z + (k) = 
(1 — <9Re£ + (u;)/<9ci;| w=ek ) _1 , and we have 



M^r 1 --^ £ ln| ek -/^|. (23) 



It exhibits a logarithmic divergence as — > //, which causes non- Fermi liquid behavior 
Eventually, N(T) is written as, 



N(T)~?f du 1- ' K 7 M UH% )7 A- (24) 



We can separate the contribution by the longitudinal gluons iVj(T) from N(T). Since 
the longitudinal gluon exchange is short-ranged by the Debye screening mass, it be- 
comes almost temperature independent, 



N,(T) = ^^fi, (25) 



with the Landau-Migdal parameter 



= ~ 8 ^ 2 / K + + «0 W - 1] , (26) 

F F 



where k = Y,f m2 D f/^kp and 



Jo(«) = ~ Z 1 : — * Jin ( - ) * ln(^ 2 ). (27) 



2 V-l 1 — M + K 
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To evaluate the transverse contribution, N t (T) = N(T) — Ni(T), we only use the 
transverse part in Eq. (|22|) : substituting Eq. ([221 into Eq. (1241) . we obtain the leading 
order contribution , 



7j-2 L 5 

C>c/ V (2fc£ 



m 



24 



kp 



-T 2 ln 



C f g 2 u F 

127T 2 



In 



0(^T 2 ), (28) 



after some manipulation. N t (T) has a term proportional to InT, which gives a sin- 
gularity at T = 0. This singularity corresponds to the logarithmic divergence of the 
Landau-Migdal parameter f( at T — 0. The chemical potential /i in Eq. fl28l) implic- 
itly includes the temperature dependence. To extract the proper temperature depen- 
dence in xm we must carefully take into account the temperature dependence of \x. 
Using the thermodynamic relation \x = —(dF/dn)\T with the free energy F = E — Ts, 
we have 1121 



M T ) = Vo- 



ir* 2kl + m 2 
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klEp 



T 1 



127T 2 



(29) 



We can see that \x includes T 2 In T term due to the dynamic screening effect for the 
transverse gluons, besides the usual T 2 term. 

As for the spin-dependent Landau-Migdal parameter, the leading-order contribu- 
tion at finite temperature comes from the transverse component / f a ; it has a loga- 
rithmic singularity at T = due to the dynamic screening effect. In this section, we 
shall see that the logarithmic divergences of iV -1 (T) and /" at T = cancel out each 
other to give a finite contribution to the magnetic susceptibility, ff is given by 



* " 27VV [1 } J (2tt)3 «9e k U ' k > 



with 



m 



An J An E s E k 



C s N- x g l M lia D t {k - q) 



(30) 



(31) 



q|= fc s 



where M na is the spin-dependent component of M u in Eq.(4), and k s = kp + 0(T 2 ) 
is defined by e^ s = \i. 

The real part of the transverse propagator is 



ReD t (k - q) 



(k - qf 



q|=fc s 



{(k-qff+dEfTTUpjmlj 



(k-q) 2 



(32) 



|q|=fc s 



: We discard here the temperature independent term of 0(g 2 ), which cannot be given only by Eq. 
However, we can recover it by taking the T — > limit later. 
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, while the imaginary part gives only a sub-leading contribution and can be discarded. 

The integral over k in Eq. (|30|) can be performed as in Eq. (21). Finally we find 
a leading-order contribution at T ^ 0, 



ft 



N~ (T) 

C f 9 2 
l2N c E s fi 



12n 2 k s E s 



InT" 1 . 



7T 2 (2k 2 



m 



-T 



InT 



-i 



0(g 2 T 2 ' 



(33) 



Compare Eq. ([331) with Eq. (27). Since E s = E F + 0(T 2 ) and k s = k F + 0(T 2 ) as 
we shall see, the InT terms cancel each other in the magnetic susceptibility ( |T9l) . 




1 2 
m(2E F + m) - -(E F + AE F m - 2m 2 )/tln- 

C f g 2 u F (2k% + ^m 2 + m 4 )^ ^ / A 



72 A; 4 ^ 2 \T 



+0(g 2 T 2 ). (34) 



In Fig. 2, we plot the magnetic susceptibility given by Eq. ( I34l) . At T=0, the 
magnetic susceptibility is positive at higher densities and the quark matter is in the 
paramagnetic phase there. At the critical density where the magnetic susceptibility 
diverges(fc^ ~ 1.6fm _1 ), there occurs a magnetic phase transition from the para- 
magnetic phase to the ferromagnetic phase and the quark matter remains in the 
ferromagnetic phase below hp. 

At T =30 MeV, there appear two critical densities at which the magnetic sus- 
ceptibility diverges. We denote these densities k F1 and k c F2 (k F1 < k c F2 ). In this 
case, k c Fl ~ 0.4fm _1 and k c F2 ~ 1.5 fm _1 . At densities below k c Fl and above k F2 , 
the magnetic susceptibility is positive, which corresponds to the paramagnetic phase, 
on the other hand, at densities between two critical densities, it becomes negative 
corresponding to the ferromagnetic phase. 

At T =50 MeV, there are still two critical densities (k Fl ~ 0.7fm _1 and k c F2 ~ 
1.3ml- 1 ), but the range between these two densities becomes narrower than at T =30 
MeV. 

At T =60 MeV, there is no longer divergence in the magnetic susceptibility and 
quark matter is in the paramagnetic phase at any density. 

We show a magnetic phase diagram of QCD on the density-temperature plane 
in Fig. 3. The four curves corresponds to the critical curves given by Eq. (l34p under 
four different assumptions: below the curves the quark matter is in the ferromagnetic 
phase, while it is in the paramagnetic phase above the critical curves. The magnetic 
transition occurs on the critical curves. 
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Figure 2: Magnetic susceptibility at finite temperature. The dotted, dashed, dash- 
dotted, and solid curves show the results at T—0, 30, 50, and 60 MeV respectively. 



For the solid curve, we have used the full expression Eq. (]34p . on the other hand, 
for the dashed, dash-dotted, and dotted curves, we have ignored the dynamic screen- 
ing^. e. the T 2 lnT term), static screening^. e. the kIuk term), and both of the two 
screenings in Eq. f !34jl respectively. 

Compare the result with the full expression flM|) with the one without the non- 
Fermi-liquid effect i. e. T 2 In T dependence. In the case without the T 2 In T term, the 
ferromagnetic phase can be sustained till over T = 60 MeV, while it can be at most 
T = 60MeV including T 2 lnT dependence. It turns out that the dynamic screening 
works against the magnetic instability and can reduce the ferromagnetic region in the 
phase diagram up to a point, but this effect is not so large. 

The dash-dotted curve is the result without the static screening or k In k term 
in Eq. (j34j) . The static screening effect works in favor of the magnetic instability to 
enlarge the ferromagnetic region. As discussed in [9], it depends on the number of 
flavors whether the static screening works for the ferromagnetism or not, which is 
peculiar to QCD. 

The maximum Curie temperature T c max is around 60MeV, which is achieved at 
kp — l.lfm -1 . Note that this is still low temperature, since T c max //cp <C 1. Thus our 
low-temperature expansion is legitimate over all points on the critical curve. One of 
the interesting phenomenological implications may be related to thermal evolution 
of magnetars; during the supernova expansions temperature rises up to several tens 
MeV, which is so that ferromagnetic phase transition may occur in the initial cooling 
stage to produce huge magnetic field. 
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Figure 3: Magnetic phase diagram in the density-temperature plane. The solid, 
dashed, dash-dotted, dotted curves show the results for the full expression Eq. (|34|) . 
the one without the T 2 In T term, without the k In k term, and without the T 2 In T 
and k\h k terms in Eq. (I34p . The open (filled) circle indicates the Curie temperature 
at kp = 1.1(1.6) fin" 1 while the squares show those when we disregard the T 2 lnT 
dependence. 

5 Outlook 

We have discussed the critical behavior of the magnetic susceptibility in the density- 
temperature plane within the Fermi liquid theory. We have found a novel non- 
Fermi-liquid behavior and phase boundary by a perturbative calculations. Some 
non-pertubative effects such as instanton effects should be taken into account at 
moderate densities. This is important not theoretically but also phenomenologically; 
more realistic estimate of the critical density or the Curie temperature is needed when 
we face phenomena in compact stars. 

There are various ideas such as amplification of the fossil field for the origin of the 
magnetic field in compact stars. So it should be very interesting if we can distinguish 
these ideas through observations. To this end we must consider not only magnetic 
evolution but also thermal evolution; if ferromagnetic state is realized, spin waves 
should be excited which affect the thermal evolution of compact stars [8]. 

This work was partially supported by the Grant-in- Aid for the Global COE Pro- 
gram "The Next Generation of Physics, Spun from Universality and Emergence" 
from the Ministry of Education, Culture, Sports, Science and Technology (MEXT) 
of Japan and the Grant-in-Aid for Scientific Research (C) (20540267). 
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